Renormalization Group Theory for a Perturbed KdV Equation 
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We show that renormahzation group(RG) theory can be used to give an analytic description of 
the evolution of a perturbed KdV equation. The equations describing the deformation of its shape 
as the effect of perturbation are RG equations. The RG approach may be simpler than inverse 
scattering theory(IST) and another approaches, because it dose not rely on any knowledge of 1ST 
and it is very concise and easy to understand. To the best of our knowledge, this is the first time 
that RG has been used in this way for the perturbed soliton dynamics. 
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The standard soliton equations(KdV, NLS, Sine- 
Gordon and so on.)are highly idealized. In more realistic 
situations, it is important to understand nonlinear evolu- 
tion equations under the influence of perturbations. 
Several different approaches to the analytic description of 
soliton dynamics in these perturbed systems are known. 
The most powerful method to deal with these cases is 
based on IST.I^"''! It turns out that for any nonlinear 
evolution equation which is reasonably close to a nonlin- 
ear evolution equation that can be exactly solved by 1ST, 
the total evolution of the scattering data can be given and 
also be expanded in a perturbation expansion. Thus one 
can determine the effects of perturbations such as small 
dissipations, relaxations, etc., on the evolution of these 
exactly solvable nonlinear evolution equations, and in 
particular, the effects on the soliton states. These equa- 
tions which describe the time dependence of the scatter- 
ing data can be solved in the " adiabatic" approximation. 
But it is inconvenient to use for one who is not familiar 
with 1ST. Another alternative way to study soliton per- 
turbations is so-called direct perturbation theory based 
on a Green's function and a rather standard 'two-timing' 
procedure. In this scheme, a basic technical ingredi- 
ent is to construct a Green's function. Then the time 
dependence of the soliton parameters and the first-order 
correction are readily available. But they stated that 
their method partially requires the knowledge of 1ST. 

In this paper we perform the RG approach to the per- 
turbed KdV equation and compare the result with an- 
other methods. We demonstrate that an analytic de- 
scription of the evolution of a soliton, including the de- 
formation of its shape and tail formation can be given, 
and the equations describing the evolution of the soli- 
ton amplitude and velocity as the effect of perturbation 



are RG equations. The apparent advantage of the RG ap- 
proach is that the starting point is simply a straightfor- 
ward naive perturbation expansion, for which very little 
a priori knowledge is required. We also see that the RG 
approach may be more efficient and concise in practice 
than another methods. 
The KdV equation 

Ut + 6uUx + Uxxx = (1) 

has a single soliton solution such as 

u^2Alsech'^z, z ^ Ao{x - 0, ^ = 4A^i + & (2) 

where Aq is the amplitude and is the initial position 
of the soliton. 

Now we consider the perturbed KdV equationl^l 



Ut + &uux + Uxxx = eR[u] 
with the initial condition 

u{x,0) = 2Alsech'^Ao{x - ^o) 



(3) 



(4) 



To proceed, we construct a e— perturbation to this equa- 
tion 



u = uo + eul + .... 
Put it into Eq.(3), we have 

Uot + GuqUox + Uqxxx = 



(5) 
(6) 



uit + GuqUix + GuoxUi + uixxx = R[uo\ (7) 



(8) 
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and so on. The initial condition now changes into 
uo{x, 0) = 2Aq sech^Ao{x — xq) 

Un{x, 0) = 0, n > 1 
In the first-order, we write 



where 



[dt- L]ui = R[uo] 



L= Guodx + 6mox + dx 



(9) 
(10) 

(11) 
(12) 



For convenience, we use z instead of x. Then Eq.(12) 
becomes 

(f d 
L= Af,[—r + (12 sec h'^z - 4) — - 24tanh^;sec/i^2l 
dz'^ dz 

(13) 

We can derive the solution to this linear partial 
differential equation (11) by the method of variables 
separation: 

ui{z,t) = ul + regular terms (14) 
so that the singular part of Ui is given by 

U'l = /lt$l(z) + (4A3/it2 + /2f)$2(z) (15) 

/ + 00 
R[uo]^i{z,k)dz (16) 
-OO 

/+00 
R[uo]Mz,k)dz (17) 
-OO 

/+00 
G{z, y; t)R[uo]dy (18) 
-OO 

1 - exp[ik{k^ + 4) Alt] 



27rfc(F + 4) 



[k^ — Aik tanh z — 4 tanh^ z] exp[iA;^;] 



4*1 (z, k) = sec h'^z 

^2{z, k) = tanh^; + zsech'^z 
Thus, the bare perturbation result is 



(23) 
(24) 
(25) 



u{z, t) = 2Al sech^z + e[ht^i{z) + {AAlf^t'' + f2t)Mz)] 
+regular terms 
z = Ao{x-0. ^ = ^Alt + ^o (26) 

which is divergent when t ^ oo as anticipatc'd. To 
deal with the divergence, we use the RG(renormalization 
group) method which has been introduced by Golden- 
fcld ot. al for global asymptotic analysis. First, to 
regularize the naive perturbation series, we introduce an 
arbitrary time r and split t as t — t + t. The quantity 
Aq and ^0. the "initial amplitude" and "initial position" 
of the soliton, can not be obtained from the knowledge 
of u{z,t) at large times, therefore Aq and are consid- 
ered to be " imobscrvablc" at large times in the same way 
that the "bare" electric charge is unobservable at long 
distances according to quantum electrodynamics. I^^l In 
another words, Aq and arc no longer constants of mo- 
tion in the presence of the perturbation. We cure this 
divergence of the bare perturbation series by introducing 
the renormalized amplitude and position 



Ar = Z^\t)Ao 



(27) 
(28) 



Giz,y;t) 

where 
^z,k) 



/ + OC 
-OO 



The multiplicative renormalization constants Zi and Z2 
are introduced to absorb the terms r. We proceed by 
assuming the expansion 



ifc(fc2 + A)Al 



$(2, fc)*(y, k)dk 
(19) 



Zi 



- [fc(fc2 + 4) -h 4i{k^ + 2) tanh z 



\/2^fc(fc2 + 4) 

—8k tanh^ z — 8i tanh^ z] exp[iA;^;] (20) 

^i{z,k) = {1 - z ta.nh z) sec h^z (21) 

^2{z,k) = tanh z sec h'^z (22) 



1 + J2^n 



(29) 



(30) 



The coefficients a„ and 6„ are chosen to eliminate the 
terms containing r order by order in e. In the first order, 
we obtain 



ai = 



fir 



(31) 
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bi = - 



.f2T 

4:A%XR 



(32) 



and we have the foUowing renormahzation perturbation 
result 



ur = 2ARsech^z + s{fi{t - t)^i{z) + 

[4AlMt-'-T-') + f2{t-r)]Mz)} + 
regular terms. 



(33) 



where Ar and ^r are now functions of r. Since r does 
not appear in the original problem, the solution should 
not depend on it. That is to say = which is 

the RG argument of Gell-Mann and LowJ^^l We get the 
following RG equation 



dAR 

dr 



1 1 /■+°° 

'4A^^'='4A^J_^ ^i^o] sec h'zdz (34) 



dr 



■4Al^ 



+ 00 



R[uo] (tanh z -\- z sec h^z)dz 
(35) 



which in this case consists of two independent equations. 
This is analogous to the RG equation in field theory. Fi- 
nally, setting T = tm Eq.(33) eliminates the secular term, 

giving the result 



u{x, t) = 2A^ sec h^z + regular term,s. 

A = Ar, z = A{x-0, ^^AAh + ^, 



(36) 



where the regular terms are the 'tail" formation!'^ ^1 and 



dA 
'dt 



'4A 



+ 00 



R[uq] sec h^zdz 



(37) 



Then we consider another example, the KdV-Burgers 
equation 



eu-r 



In the same way, we obtain 



A = 



Ao 



16^/12 



C=glog(l + i^..A^) 



(41) 



(42) 



(43) 



which is also the same as that derived by IST.I^^^ 

In summary, we have demonstrated that a perturbed 
KdV equation can be solved by RG method, with some 
attendant technical advantages compared with other 
methods. The present approach can easily be general- 
ized also to multiple-soliton state. To avoid unnecessary 
complications we expound our theory using as example 
the KdV equation. It is, however, clear from what fol- 
lows that this does not restrict the general nature of the 
method for another perturbed evolution equations. 
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!^=4A'^-^£—— R[uo\{ta.nhz + zsech'^z)dz 

(38) 

The two important equations which determine how the 
soliton shape and position are affected by the perturba- 
tion are also consistent with those derived by IST.['^~^i 

As an example, we consider the damping KdV equa- 
tion in which R[v] = —u. The time dependence of the 
soliton parameters can be easily obtained from Eq.(37) 
and Eq.(38). Thus we have 

A = AoeM-'^) (39) 

e=^^^[l-exp(-^)] (40) 
which is just the same as that obtained by IST.!^''! 
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